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Abstract 

The effect on the transport properties of fractures of a relative shear displacement 
u of rough walls with complementary self-affine surfaces has been studied experi- 
mentally and numerically. The shear displacement u induces an anisotropy of the 
aperture field with a correlation length scaling as u and significantly larger in the 
direction perpendicular to u. This reflects the appearance of long range channels 
perpendicular to u resulting in a higher effective permeability for flow in the direc- 
tion perpendicular to the shear. Miscible displacements fronts in such fractures are 
observed experimentally to display a self affine geometry of characteristic exponent 
directly related to that of the rough wall surfaces. A simple model based on the 
channelization of the aperture field allows to reproduces the front geometry when 
the mean flow is parallel to the channels created by the shear displacement. 
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1 Introduction 



The geothermal reservoir of Soultz-sous-Forets, like most geological systems, 
contains fracture networks providing preferential paths for fluid flow and iden- 
tifying these paths is crucial in geothermal applications. More precisely, the 
number of hydraulically connected fractures and their spatial distribution are 
key factors in determining the rate at which hot water can be produced or 
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fluids can be re-injected into the system ( Center et al 1997 : Dezaves et al 
200J). Transport in these systems can be stimulated by the massive injection 
of water at high flow rates which is commonly used to increase the overall 
transmissivity : both borehole and microseismic observations have demon- 
strated that t he relative shearing of the fracture wa lls is a key factor in this 
enhancement ( Center et al 1997 : Evans et al 2005[ ). Another important is- 
sue is the influence of the shear on the transport of suspended particles and 
dissolved species. 



The present work analyzes the influence of shear displacements of the fracture 
walls on the pore space structure and, therefore, on the transport properties. 
The case of fractures with a self-affine geometry is speciflcally investigated 
this statistical mo del describe s inde ed successfully a broad variety of frac- 



tured materials (see Bouchau'dl ( 2003| ) for a review). This model which will 



be described more precisely in section 3 accounts also well for the long range 
correlations in the topography of such surfaces. 

An important factor in these processes is the anisotropy of the permeability 
fleld induced by the shear displacement of the fracture walls : the permeability 
is highest fo r flow perpend i cular to the shear disp lacement vector u and low- 



est along it (ICentier et al 



^ , , 119971: Yeo et al. . 1998| ). This anisotropy has been 

suggested by Auradou et all ( 2005 ) to be associated to the appearance, in the 
aperture fleld, of ridges perpendicular to u. These ridges partly inhibit flow 
parallel to u while they create long preferential flow channels connected in 
series in the perpendicular direction : this leads to model the pore space as 
long parallel channels extending over the length of the fracture. In the present 
work, the validity of this simple description of the void space will be conflrmed 
by a detailed study of the geostatistical characteristics of the aperture fleld. 
It will be shown that both the anisotropy of the aperture fleld and the pres- 
ence of large scale surface structures spanning the whole surface are needed 
to account for the anisotropy. 



These analysis will then be extended to the spreading of a tracer advected by 
the flow. This process is largely determined by the structure of the flow fleld 
which w ould, at flrst, seem t o be very complex. However, both experinienta. 1 
studies ( Auradou et a and numencal simulations JtoiZSIi» 

suggest that, for shear displacements m small compared to the fracture size and 
for flow perpendicular to u, the flow lines are mostly parallel to the mean flow 
the velocity fleld remains correlated over distances of the order of the fracture 
length. Finally, the correlation of the velocities of tracer particles advected in 
the fracture will be shown to have also self-affine characteristics, and a model 
accounting for the geometry of miscible displacement fronts in such fractures 
will be presented. 
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Fig. 1. Schematic view of the experimental fracture models. The fracture walls are 
complementary surfaces with the same roughness geometry and self-affine exponent 
C ~ 0.8 as natural granite fractures. The model walls are shifted vertically to deter- 
mine the mean aperture and horizontally to model relative shear displacements. 

2 Experimental setup and numerical models 



Both in the experiments and the numerical simulations, the fractures are mod- 
eled as the gap between two complementary rough surfaces, as shown in Fig. 1. 
Experimentally, these surfaces are transparent casts of granite or fracture sur- 
faces that retain their self affine geometry. The experimental setup allows us to 
control both the normal distance and the lateral displacement b etween the two 
facing wall surfaces (a detailed description has been given bv lAuradou et 



( 2001 . 2OO5I ) ). The normal distances investigated range from 500/xm to 1 mm, 



i.e. of the same order of mag nitude as the mean aperture of fractures in the 
Soultz Geothermal reservoir (|Saussel . l200i . The lateral displacement mimics 



a relative shear displacement of the walls and induces spatial variations of 
the aperture. The magnitude of the shear displacements ranges from zero (no 
shear configuration) to a few millimeters when the lateral displacement brings 
the opposite surfaces in contact. The shear displac ement ampli t udes are of the 
same order of magnitude than those estimated by lEvans et a after an 



hydraulic stimulation. Th e result i ng ap erture fluctuations are also similar to 
those reported by Sausse ( Saussd . 2002) for fractures in the Soultz reservoir. 



The fracture model is fitted inside a parallelepipedic basin : both are initially 
carefully saturated with a transparent solution so that no air bubble gets 
trapped inside the setup. The same solution, but dyed, is then injected through 
a small borehole in the middle of the upper epoxy cast and pictures of the dye 
invasion front are obtained at regular time intervals. Figure 2 displays a typical 
invasion front geometry : its larger elongation in the direction perpendicular to 
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Fig. 2. a) Experimental pictures of the radial injection of a dyed fluid into a fracture 
saturated by a miscible, transparent one (picture size : 20 cm x 20 cm. b) Detail of 
the invasion front. 

the shear reflects the influe nce of the permeabiht y anisotropy discussed above. 
Moreover, as discussed bv lAuradou et a the ratio of the elongations 

of the front parallel and perpendicular to the shear is directly related to the 
ratio of the permeabilities along these two same directions. 



As mentioned above, aperture variations due to the lateral displacement of the 
fracture walls strongly perturb the flow fleld. The numerical simulation of the 
above experiments requires therefore an accurate 3D description of the frac- 
ture geometry and of the velocity fleld : the lattice Boltzmann technique has 
been selected for that purpose becau se it is well adapted to complex geome- 
tries (|Drazer and Konlikl . Efflnl 120021 ) . A large cubic lattice (1024 x 1024 x 20 
sites) is used to reproduce a model fracture of mean aperture 1 mm and area 
51.2 X 51.2 mm (one fourth of the experimental area). Using these numerical 
simulations, the variability of the results from one fracture to another has been 
investigated and the permeabilities for flow both parallel and perpendicular 
to the shear direction have been determined (note that only their ratio can be 
measured experimentally). 



Figure 3 compares aperture flelds obtained in numerical simulations and cor- 
responding to relative displacements of the same complementary surfaces in 
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Fig. 3. Numerical aperture fields (51.2 mm x 51.2 mm) obtained for a same mag- 
nitude of the shear displacement u {u = 0.8 mm) along two different perpendicu- 
lar directions (arrows). Upper images - aperture fields coded in grey levels. Lower 
images : binarized images corresponding to a 1 mm threshold equal to the mean 
aperture (lower apertures appear in white). Left : S = 0.16, K^/Kq = 1.05 and 
K\\/Kq = 0.822 {S = ratio of the variance of the aperture fluctuations by the 
mean aperture, K\\ and : permeability values for flow parallel and perpendicu- 
lar to the shear displacement n, Kq = permeability for u = 0.). Right : S = 0.21, 
K^/Ko = 1.00 and K\\/Kq = 0.930. 

two directions at right angle. Permeability values computed in these two cases 
for constant pressure gradients parallel and perpendicular to the shear are 
listed in the caption. One observes the appearance of near continuous streaks 
perpendicular to the shear displacement and marking preferential flow paths. 
The appearance of these paths is reflected in the difference between the per- 
meabilities i^x and K\\ for flow respectively parallel and perpendicular to 
them. More precisely, the permeability K± for flow perpendicular to the dis- 
placem ent u increas es with u while K\\ decreases. These variations are ob- 
served (jAuradou et al . .20051) to be related to the permeability Kq for m = 



5 



and to the ratio S of the variance of the aperture fluctuations, cTq by the mean 
aperture, Oq through the approximate equations : 



^ = 1-6BS' + 0{S'), (1) 



The two constants, A and B, depend on the specific statistical properties of 
each individual fracture surfaces but satisfy the relations: A < 1 and B < 1 . 



As stated above, the binarized aperture fields shown in Fig. 3 displays ridges 
perpendicular to the dispacement u and extending over a large part of the 
fracture length parallel to the flow. This observation has been used as the 
basic ingredient of a simple model in which the aperture field is replaced by 
a set of parallel channels, each of constant aperture and perpendicular to the 
shear. This model captures the linear dependence of the permeabilities with 
respect to 5^ and giv es a good estimate of the respective upper and lower 
bounds of K i and K h ( Zimmerman et al 1991 : Zimmerman and Bodvarsson . 
Il996l : lAuradou et al I2nn5[ ). These two limits correspond to A = 1 and B = 1 
in Eq. (1). This simple model takes advantage of the fact that the correlation 
lengths in the aperture field are highly anisotropic and much smaller parallel 
to u than perpendicular to it. 



In addition to its global anisotropy, the experimental pictures displayed in Fig- 
ure 2 show that the front between the transparent and dyed fluids is locally 
quite tortuous. More precisely, its geometry has been found experimentally 
to be self-affine with a c haracteristic exponent close to that of the roughness 
of the fracture surfaces (lAuradou et al 2001 ). This result has been verified 
numerically bv lDrazer et a ll (I2OO4I) for two values C = 0.8 and 0.5 of the char- 
acteristic exponent ( (Feder|"^8^ of the fracture wall surfaces (note that C 
has the same "universal" value 0.8 for the basalt and granite surfaces mod- 
eled in the present work and its influence could not therefore be investigated 
experimentally) . 



The spatial correlation of the aperture and velocity fields in these fractures 
will now be analyzed more quantitatively using standard geostatistical tools. 
A model based on these caracteristics will then be discussed in section 4 and 
shown to capture the key features of miscible displacement processes studied 
experimentally and numerically. 
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3 Spatial correlation of the aperture field 



The following interpretations are realized on 2D maps with 770 x 760 pixels 
measured experimentally on the epoxy surfaces by a mechanical profilometer; 
the spacing between pixels is 250 /xm in both directions of the mean plane of 
the fracture and the vertical resolution on the local height h{x, y) is 10 /im. 
The statistical properties of the surface roughness have determined from these 
maps and, as stated above, can be characterized in terms of a self-affine geom- 
etry. The main property of a self-affine surface is that it remains statistically 
unchanged under the transformation : 



X — > Ax, 
h{x,y)->\^h{x,y). 

Where (x, y) are the coordinates in the mean surface plane, y) is the 
local height of t he surface ari d C, is refered to as the Hurst, roughness, or self- 
affine exponent ( Fedei . 1988[). For the experimental model fractures, one has 
C = 0.75 ± 0.05 (jAuradou et a 1, l2005f), in agreement with published values 
for granite surfaces (Bouchaud . .2003^). In the simulations, the self-affine sur- 
faces are generated numerically (see 
description of the method). 



Drazer and Koplik (|2002f ) for a detailed 



In the present work, we consider fractures made of complementary walls that 
are separated both vertically in order to open the fracture and laterally in 
order to mimic shear displacements (if any). Under such conditions, if oq is 
the mean aperture and u the relative shear displacement in the plane of the 
fracture, the local aperture at a point f satisfies : 



a(r) = h{r) — h{r + u) + oq. 



(2) 



This equation allows to determine the aperture field a(r) for a given shift u 
from the experimental surface profile map h{r). The correlations in the aper- 
ture field are characterized by the following correlation function, also called 
semivariance (|Kitanidij . lT997[) : 



7(5) = ((a(f)-a(f+5))2), 



(3) 



which measures the spatial correlation of the aperture field between two points 
separated by the 2D vector 5. 

Combining Eqs. 2 and 3, one sees that the value of 7 is determined by the 
correlation between the heights of the surface at four points located at the 
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Fig. 4. Semivariograms of aperture fields corresponding to shift amplitudes u equal 
to : (0.25 mm, o), (0.5 mm,n), (0.75 mm,0), and (1 mm, Dotted lines : 7|| 
(correlation along the direction of u). Solid lines : 7^ (correlations in the direc- 
tion perpendicular to u). Horizontal dashed lines correspond to twice the aperture 
variance for the same u values. 



corners of a parallelogram of sides u and 5. The orientations of the lag 5 
respectively parallel and perpendicular to the shear displacement u are of 
special interest (in the following, they are respectively referred to as || and _L). 



Figure 4 displays variations with 5 of the semivariances 7^ and 7|| for aperture 
fields computed from a same experimental surface profile map h{r) but for 
different shear displacements u. For small values of the lag S, the apertures at 
r and f+6 are very similar (high spatial correlation at short distances) and the 
value of 7 is low. As S increases, so does 7 because the correlation decreases. As 
6 increases further to a value higher than the correlation length of the aperture 
field, the apertures at r and f+6 become uncorrelated; then, from Eq. (3), 
7_L and 7|| should reach a same constant saturation value equal to twice the 
aperture variance tXa = ((a(x, y) —aoY). The values of aa corresponding to the 
different shear amplitudes u are plotted as horizontal dashed lines on Fig. 4 
: they represent ind eed the limit of both 7 and 7|| at high 6 values. Also, as 
previously shown bv lAuradou et a (Jq increases with the magnitude u 

of the shear displacement : in the range of u values used in Fig. 4, one verifies 
experimentally that cTq = Au'^'' with a best fit corresponding to ( = 0.75 and 
A = 137000 mm2(i-?). 
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Fig. 5. Semivariograms of aperture fields normalized by Au^"^ {(^ = 0.75 and 
A = 137000 mm^(^~^ as function of the normalized lag 6/u. Shear amplitude u 
values : (0.25 mm, o), (0.5 mm,\3), (0.75 mm,0), and (1 mm, *). Dotted lines : 7|| 
(correlation along the direction of u). Solid lines : 7^ (correlations in the direction 
perpendicular to u). Note that the horizontal range is narrower than in Fig. 4 so 
that the overshoot effect is not visible. 

Even though the limits of 7^ and 7|| are the same at large 5 values for a given 
shear displacement u, the way in which this saturation is reached differs very 
much for the two semivariances due to the anisotropy of the aperture field. The 
semivariance 7_l for 6 perpendicular to u never exceeds its saturation value 
and reaches it in an overdamped way whereas 7|| {5 parallel to u) displays an 
overshoot reflecting a local anticorrelation. These different behaviors result 
directly from the large scale anisotropic structures observed above on Fig. 3. 
Let us consider the aperture at two points separated by a small distance 6 
perpendicular to the shear displacement (and therefore parallel to the length 
of the structures) : they will most hkely be both either above or below the 
average aperture if they are inside a same large scale structure. Then, the 
corresponding value of aperture semivariance will be less than (Tq. On the other 
hand, two points with a separation 5 parallel to u have a high probability to 
correspond to different structures of the binarized picture : some amount of 
anti correlation can therefore be expected. 

Let us now investigate the parameters determining the variation of 7 in the 
transition regime. For that purpose, we have seeked to collapse the different 
curves of Fig.4 onto master curves. In Fig. 5, 7^ and 7|| are normalized by 
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Fig. 6. Variations of (7||/7|| — 1) (top) and (7_l/7± — 1) (bottom) normalized by 
u'^ (C = 0.75) as a function of 6. Shear amplitudes u : (0.25 mm, o), (0.5 mm,n), 
(0.75mm,O), and (1 mm, *). 

Au'^'' which represents their limit at large va lues of 5 while, for the ho rizontal 
scale, 6/uisa natural reduced coordinate (see Plouraboue et al. ( 1995[ )). These 
two normalizations allow indeed to collapse in Fig. 5 the different variations 
into two independent master curves, corresponding to the two orientations of 
6 with respect to u. This demonstrates, that, in all cases, the characteristic 
distance over which the aperture field gets decorrelated is proportional to the 
shift u : in addition, the ratio of this distance by u is significantly larger for 
6 perpendicular to the shear than parallel to it. While this latter result will 
remain valid, the actual values of both distances may vary significantly from 
one sample to another. 



Fig. 5 also shows that, for small lag values 6/u -C 1), the semivariograms 
increase linearly with 6/u {j/u^ oc 6/u). This indicates that the random field 
is not as smooth as in the case of a Gaussian covariance function for which 
7 oc 6"^. 

At large distances 5 ^ m, the semivariance 7 still displays large fluctuations 
with excursions above and below the asymptotic value cXa : these are visible in 
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Fig. 6 that represents the normahzed deviations (7— 7)/7 of 7 from its average 
value 7 over all distances 6 > 15 mm. It is observed that the fluctuations 
coincide provided 6 is used as the horizontal scale (and not the normalized 
lag as in Fig. 5). This suggests that these variations reflect the large scale 
underlying structures displayed in Fig. 3 and that their location is independent 
on the magnitude u of the shear displacement (these fluctuations will however 
be very variable from a sample to an other one). The variation of the amplitude 
of the fluctuations with u is however non trivial : the normalized deviations 
(7 ~ 7)/7 had to be divided by u'^ to obtain the collapse of Fig. 6. This 
indicates that there is some influence of the local self-affine structure on these 
fluctuations; if they were only due to the large scale channels, the normalized 
deviations would be proportional to u. 



4 Correlation between fracture aperture fields and fluid displace- 
ment fronts 



This section is devoted to the geometry of the front between two miscible flu- 
ids (or equivalently to the dispersion of a tracer) in a single fracture. Several 
mechanisms contribute to the spreading of the front (or of the tracer). These 
include molecular diffusion, Taylor dispersion related to the velocity profile 
in the fracture gap and geometrical dispersion due to velocity variatiori s be- 
tween flow li r ies in the mean fracture plane fsee lDrazer and Koplik (HooJ and 



Drazer et all ( 2004 ) for a discussion of the contribution of the above mentioned 



mechanisms to tracer dispersion). 

Here, we shall take only into account the last mechanism : this allows us to 
assume a two-dimensional flow field v{x,y) obtained by averaging the three- 
dimensional velocities over the gap of the fracture with v{x, y) = {v{x, y, z))^. 
The pressure gradient inducing the flow and the resulting global mean velocity 
are parallel to the axis x and the front (or the tracer) is initially located at 
the inlet of the fracture on a line parallel to the axis y. In the following, we 
analyze the development of the front (or the tracer line) with time, assuming 
that its various points move at the local flow velocity, and we present a simple 
analytical model accounting for its geometry. 

This model is based on the results discussed above, namely that the aperture 
field is structured into channels perpendicular to the shear displacement u. For 
a mean flow parallel to these channels, the tracer particles follow paths that are 
only weakly tortuous; also, the variations of their velocity along these paths are 
small compared to the velocity contrasts between the different channels. Under 
these assumptions, the velocity of a particle located at a distance y measured 
perpendicular to the mean velocity satisfies : v{x,y) ~ v{y)nx, where fix is 
a unit vector parallel to the mean flow. Note also that there are no contact 
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Fig. 7. Compared normalized front profiles (circles) and normalized apertures 

(o^Hy)) lines) as a function of the distance y. across the flow. Top and middle 

graphs : flow normal to shear displacement u with u = 0.2 mm and u = 0.6 mm. 
Bottom graph : Flow parallel to u for u = 0.6 mm (the y scale of the solid line has 
been amplified by a factor 10). For all graphs, the fracture aperture is oq = 1 mm 
and its size is 51 mm x 51 mm. 

points between the walls of the fractures we have used : this avoids to take 
into account the large tortuosity of the flow lines in their vicinity. For each 
channel the velocity v{y) is related to the pressure gradient dP/ dx by Darcy's 
equation : 



a^jy) dP 
12yU dx 



(4) 



in which a{y) is the equivalent (or hydraulic) aperture of the channel. Then, 
particles starting at different distances y at the inlet move at different veloc- 
ities, v{y), and their distance x from the inlet at time t after the injection 
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satisfies x{t) = v{y)t; the mean front position moves tlien at tlie average ve- 
locity witli x{t) = {v{y))yt. Combining tlie previous relation with Eq. (4) we 
obtain: 

x{y,t) ^ a^{y) 

x{t) {a\y))y- ^ ^ 



Finally, previous studies have shown that, for relatively small aperture fluctua- 
tions, the hydraulic aperture is we l l appr oximated by the geometrical aperture 
f Zimmerman et al 1991 : Brown . 1987 ) : this suggests that a{y) can be ap- 



proximated by the average of the local apertures along the direction x with 

a{y) = {a{x,y))^. 

The validity of this assumption has been tested by comparing the shape of the 
front x{y, t) determined from simulations using the lattice Boltzmann method 
to the variations of the square a^iy) of the effective aperture. For that purpose, 
the normalized displacements x{y, t)/x{t) at a time t parallel to the mean flow 
of particles released initially at x = are plotted in Fig. 7 as a function of 
the distance y together with the normalized square of the aperture averaged 
parallel to the flow : a^{y)/ {a'^{y))y. When the mean flow is perpendicular to 
the shear displacement u, the variations with y of the normalized position of 
the particles and of the normalized square aperture are very similar both for 
u = 0.2mm and u = 0.6mm (top and middle curves). The small differences 
may be due in part to the effect of viscous entrainment between adjacent flow 
channels due to viscous diffusion. On the contrary, when the mean flow is 
parallel to u (bottom curves), there is very little correlation between the two 
curves and the amplitudes of the variations are very different (note that the 
curve corresponding to a?{y)/ {a?{y))y has been magnified by a factor of 10). 

These strong differences reflect and confirm the anisotropy of the aperture 
field induced by the shear displacement u : the assumptions underlying the 
above model and, particularly, Eq. 5 are indeed only valid if the flow field 
can be described as parallel channels with a length of the order of that of 
the fracture. For a mean flow perpendicular to m, this model represents a 
good approximation : the aperture of the flow channels displays long range 
correlations and particles are convected along each channel with a weakly 
fluctuating velocity. Note that, in this case, the global width of the front 
parallel to x increases linearly with time while, for dispersive processes it 
increases from the square root of time. This difference results from the fact 
that the correlation length of the velocity of fluid particles is of the order of 
the total path length through the fracture; a dispersive regime can only be 
reached if the correlation length is much smaller. For flow parallel to m, the 
flow lines are more tortuous and velocity fluctuations are enhanced : the model 
does not account any more for the geometry of the front and its width parallel 
to the mean flow is in addition reduced. 
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Fig. 8. Log-log plot of the power spectra P{f) of the deviations of the front dis- 
placement x{y,t) (dashed line) and of the squared effective aperture a^(y) (con- 
tinuous line) as a function of the spatial frequency / = L/x (shear displacement 
u = 0.2 mm, mean aperture oq = 1 mm, fracture size L = 51.2 mm). The curves 
are shifted vertically relative to each other for convenience. 

Further informations are obtained by performing Fourier transforms of the 
deviations of both the displacement x{y,t) of the front and the square a^{y) 
of the estimated effective aperture from their mean values (the mean flow is 
perpendicular to the shear displacement u). The resulting power spectra are 
displayed in Fig. 8 and, as could be expected from the qualitative analysis of 
the curves of Fig. 7, both spectra arc similar up to frequencies / of the order 
of / = L/u. In fact, for wavelengths smaller that the shear displacement u, 
only a'^{y) displays self-afRne properties. 



5 Discussion and Conclusions 

The experiments and numerical simulations realized in the framework of the 
present program have revealed new features of flow and transport in rough 
fractures which may be of importance in the modelization of fractured field 
reservoirs, such as those encountered at the Soultz site. We have specifically 
studied the influence of a relative shear displacement of complementary wall 
surfaces on the aperture fleld as well as on flow and transport between them. 
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This is particularly relevant to geothermal reservoirs where fluid injections 
frequently induce relative displacements of the fracture walls. In agreement 
with many experimental results, the rough fracture walls have been assumed 
to display a self-affine geometry. 

First, we have shown that the relative shear displacements -u induce an anisotropy 
of the flow field which reduces permeability for flow parallel to u and enhances 
it in the perpendicular direction; this may account in part for the permeabil- 
ity enhancement in some sheared zones of the reservoirs. These variations are 
well predicted by assuming the appearance in the aperture field of large ridges 
perpendicular to the shear and extending over the full length of the fracture 
: these ridges act as parallel channels enhancing flow parallel to them and 
reducing transverse velocity fluctuations. This anisotropy has been confirmed 
quantitatively by computing semivariograms of the aperture fields : these dis- 
play a correlation length proportional to the shear displacement u for a given 
orientation but significantly larger perpendicular to u than parallel to it. 

A second important result is the strong correlation between the roughness of 
the fracture walls and the geometry of the displacement fronts, particularly for 
fiow perpendicular to u. In particular, the displacement front has a self-affine 
geometry characterized by the same exponent as that characterizing the frac- 
ture walls. Also, for a mean fiow perpendicular to u, the front width is found 
to increase linearly with distance : this refiects a channelization parallel the 
mean fiow, resultinb in weak velocity fiuctuations along the paths of the fiuid 
particles. Another important consequence is the fact that the front geometry 
can be estimated accurately in this case from the variations of the square of 
the aperture averaged parallel to the mean fiow. 

These results raise a number of questions that will have to be considered in 
future studies. First, one may expect the spatial correlations of the velocity 
field to decay eventually, leading to a normal Gaussian dispersion process at 
long distances. This has not been observed in our experimental or numerical 
systems and wc have currently no estimate of the corresponding correlation 
length. Also, the channeling effect might result in a motion of some fiuid 
particles at a velocity substantially faster than the average fiow velocity, at 
least in transient regimes of duration shorter than the time necessary to sample 
the fiow variations inside the fracture. This effect may have a strong influence 
on the thermal exchange between the rock and the fluid, or on the transport 
of species in dissolution-deposition processes encountered, for instance, on the 
Soultz site. In the case of a channelized fiow, for instance, the variability of 
the residence times in the various fiow channels may result in temperature 
inhomogeneities in the outfiowing fiuid. 
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